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SUMMARY 

This design guide describes a method for predicting the random response cf flat 
and curved plates which is based on theoretical analyses and experimental results. 
The plate curvature can be due to postbuckling in-plane mechanical or thermal 
stresses. Based on a single mode formula, r.m.s. values of the strain response to 
broadband excitation are evaluated for different static buckled configurations using 
the equivalent linearization technique. The effects on the overall strain response 
due to instability motion of snap-through are included. Panel parameters include 
clamped and simply-supported boundaries, aspect ratio, thickness and length. 
Analytical results are compared with experimental results from tests with 12 in. x 
15 in. aluminum plates under thermal loading in a progressive wave facility. 
Comparisons are also made with results from tests with a 2 in. x 15 in. x 0.032 in. 
aluminum beam under base mechanical excitation. The comparisons help to assess 
the accuracy of the theory and the conditions under which deviations from the theory 
due to effects of imperfection and higher modes are significant. 

INTRODUCTION 

Although much work has been done recently on developing analytical methods 
for predicting nonlinear random response of structures (ref. 1-4), most of these 
procedures are too tedious for practical application and have not been 
experimentally validated. Therefore, it is necessary to rely on empirical formulae 
(refs. 5-6) and experience to design aircraft structures against sonic fatigue. 
Furthermore, the estimation of nonlinear responses in plate structures (especially 
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curved plates) is so uncertain that designers often choose to avoid nonlinear 
problems by making thicker and heavier structures. The resulting increase in weight 
may be unacceptable in many future classes of high performance aircraft. On the 
other hand, it should be noted that for flat and curved plates, nonlinear response can 
result in hardening-spring effects (increase of stiffness with amplitude of deformation) 
and thus offer weight saving if the dynamic deflection can be controlled. 

Some of the reasons for the lack of validation of prediction methods for 
nonlinear random vibration due to acoustic loads are the difficulties in generating 
large enough forces with desired spectra and difficulties in maintaining consistent 
rigid boundaries of plate specimens without premature fatigue problems. The 
availability of high intensity acoustic test facilities is limited and most acoustic tests 
have been performed on built up structures (ref. 5-9). Although these tests have 
recorded large amounts of structural response and fatigue life data, no vigorous 
attempts have been made to compare these data with theoretical predictions. There 
is also a lack of experiments on nonlinear response of curved plates. Design 
guidelines have been suggested but these are only applicable to structures which 
are similar to the test specimens and in general neglect effects of nonlinear 
response. 

The status of nonlinear random response research is that there is a significant 
amount of purely theoretical analyses and some practical experimental work but little 
or no interaction between the two. To bridge the this gap, simplifications in the 
procedures on both sides are required. The compromising solution and the 
approach taken in this reported activity is to derive simple and approximate 
analytical prediction techniques and improve these through comparison with 
experiments on simple structures, such as, beams and single bay plate specimens. 

A Rayleigh-Ritz analysis using one mode for representing the transverse 
displacement of a plate was developed to predict the buckling displacement due to 
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in-plane stresses and the resulting equation of motion. The nonlinear acoustic 
response including the mean and r.m.s. strain response are predicted by the 
equivalent linearization technique. Previous work (ref. 10) using this technique has 
been shown to give good correlation with experiment. The prediction procedure is 
capable of handling hardening and softening spring effects and snap-through motion 
from one static position to another. Snap-through motion has significant effects on a 
buckled plate and thus it is necessary to evaluate the ability of the sound pressure 
level to initiate snap-through motion and to evaluate overall r.m.s. strain response 
with frequent snap-through. Design charts have been compiled and are presented 
in this report. Comparison of analytical and experimental results have been made 
with data from acoustic tests on a thermally buckled plate and from a mechanical test 
with a buckled beam. 
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SYMBOLS AND ABBREVIATIONS 

width of the plate, in. 
length of the plate, in. 
frequency, Hz 

fundamental resonance frequency 

acoustic pressure, lb/in 2 

linear modal stiffness of the flat plate, lb./in 2 

nonlinear modal force coefficient, lb./in 3 

nondimensional force parameter, = P/KQn 

, a »b 

externally applied modal force P= J J h <f>(x, y)dy dx 

o o 

P = KQ + KnQ 3 for a flat plate 


qo 

Q 


nondimensional displacement coefficient = Q/Qn 

value of q for static buckling 

modal displacement (see definition w below) 
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Qn 

R 

s 

Shh(f) 

Spp(f) 

Sp P (f 0 )L or S L 


. S s 
S c 

t 

u 

Uc 

w(x,y) 


value of Q for which linear restoring equals the nonlinear restoring 
force i.e. KQn = Kn Qn 3 - Note also Q=Qn when R=1 

u/u c - 1 (the buckling parameter) 

aspect ratio = b/a 

spectral density of h at frequency f, dB/Hz 

spectral density of p at frequency f, Hz -1 

value of Spp(f) when r.m.s. response of q is 1, assuming no 
nonlinear restoring force S pp (f 0 ) L = 

value of Shh(f) when snap-through motion starts, dB/Hz 

value of Shh(f) when continuous snap-through motion starts, dB/Hz 

thickness of the plate, in. 

inplane edge shortening displacement, in. 

value of u at which buckling starts, in. 

transverse displacement of the plate at coordinate (x,y), in. 
w(x,y) = Q <Kx,y) 


x.y 


s pp< ( > 

a=S pp«o> L 

Q 

c 

p 

Y 

e 

£0 

£N 

<>(x.y) 


coordinate system with x-axis and y-axis parallel to shorter side, a, 
and longer side, b, respectively 

nondimensional excitation parameter 


linear circular frequency of the flat plate, radians 
modal damping coefficient 
density of the plate material, lbm/in 3 
in-plane stiffness parameter 

~2 

surface strain in microstrain due to bending = -i -2J" 

2 dx 2 


value oT e for static buckling 

value of e corresponds to Qn 

shape function of the fundamental mode, 


0 


sin^r- sin-^- for simply supported boundaries, 

a b ' 

T(x/a) ^(y/b) for clamped boundaries 
where H'(z) = 5/8 [coshdz - cosdz - 0.98 (sinhdz - sindz)], d = 4.73 
Note that maximum <J)(x,y) is 1 

ANALYTICAL MODEL FOR SNAP-THROUGH MOTION 

Most previous work (ref. 1 1, 12) on theoretical prediction of snap-through under 
random excitation was concerned only with the stability boundaries but not with the 
magnitude of the nonlinear response. Also, the methods developed are lengthy and 
not easily applied to practical case of plates under acoustic excitation. 

From reference 13, it is found that under low level random excitation a plate 
exhibits snap-through motion occasionally but the oscillation is generally around the 
initial static position. This is defined as the intermittent snap-through. Continuous 
snap-through motion is defined for higher levels of random excitation when snap 
through motion constitutes more than 80 percent of the response. This usually 
occurs when the r.m.s. value of the dynamic displacement is equal to the initial static 
displacement. The general characteristics of the three regions of snap-through 
motions are summarized in Table 1. 

The division of snap-through motion into three regions helps to understand the 
instability motions of buckled plates. The boundaries between them, however, are 
sometimes difficult to define, especially for shallow buckled plates, due to the small 
gap between the two static positions, or for buckled plates which have only one static 
position. 

Snap-through boundaries were estimated empirically in reference 14 from a set 
of experimental results but no general formulae were derived. A formula for buckled 
plate behavior and an analysis model for the random response based on the above 
three types of snap-through motion was established in reference 1.0 and will be 
discussed next. 
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General Equations for Nonlinear Behavior of Plates 

Assuming an appropriate shape function for the transverse displacement of a 
plate and using the Rayleigh-Ritz procedure of formulation, the nonlinear equation 
for equilibrium in the transverse direction is obtained. The detailed steps in deriving 
the following equations are found in reference 10. 

For a plate under uniaxial compression with uniform edge displacement, the 
nondimensional relationship between modal displacement, q, and modal force, p, for 
the buckling mode is given as: 
for static equilibrium 

q 3 - Rq = p (la) 

and for dynamic motion 

q/Q 2 + 2£q/Q + q 3 - Rq = p (15) 

where R is the buckling parameter given by u/u c -1, u is the in-plane edge 
displacement, u c is the critical value for initiation of buckling, Q is the frequency of the 
plate and £ is the modal damping coefficient. The above equations are applicable to 
any plate of any boundary condition and aspect ratio if there is no coupling between 
the various modes. A plot of p versus q for various values of R from (la) is shown in 
figure 1. From figure 1 regions of hardening and softening spring behavior and also 
regions of negative stiffness, (e.g. between A and B) can be found. Notice that for 
R=1, dynamic motion starting from C will pass through A and B and end up at C. This 
motion through the instability region is called snap-through. A and B are found to be 
at (0.577, -0.385), (-0.577, 0.385) respectively by locating the points where dP/dq = 


0 . 
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Static equilibrium positions q 0 for positive values of the buckling parameter R 
are found by putting p=0 in equation (la) and correspond to the points where the 
curve crosses the q axis in figure 1 , qo=^^ - ' and zero. The last value, zero, 
is an unstable position since the stiffness is negative. Equations (la and 1b) can be 
written in terms of the postbuckling displacement q 0 as 


„3 _2 _ _ 

q -q 0 q = p 

(2a) 

q/n 2 + 2Cq/<T + q 3 -q 3 q = p 

(2b) 


Equations for Predicting Nonlinear Random Responses 

Since the mean values of snap-through motion changes with increasing 
excitation level, the technique of equivalent linearization is used to evaluate the 
mean value as well as the mean square value of the random response. The 
procedure is similar to that used for a problem of random vibration of cylindrical 
plates (ref. 15). 

Substituting q = q+ Aq where qis the mean value and Aq is the oscillating 
component and noting that Aq=0, it can be shown that 

q 3 = (q + Aq) 3 = q 3 + 3 g 2 Aq + 3q(Aq) 2 + (Aq) 3 

. = q 3 + 3q(Aq) 2 _ < 3 > 

Since P = 0, and using the value of q 3 from equation (3) the mean of the equation 

of motion (2a) is given by 

q 3 + 3 q(Aq) 2 - q 2 q = 0 
The solutions are 

q 2 = q 2 - 3(Aq) 2 andq = 0 

Using equivalent linearization, the dynamic equation (2b) can be written as 


( 4 ) 
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Aq + 2Q£Aq + Q 2 (3q 2 - q 2 )Aq - £l 2 p 


( 5 ) 


Substituting q 2 = q 2 + (Aq) 2 and q 2 = q 2 - 3Aq) 2 from equation (4) into (5), the 


dynamic equation becomes, 


Aq + 2Q£Aq + Q 2 [2q 2 - 6(Aqf] Aq = Q 2 p 


( 6 ) 


From equation (6), the random response is given by (See Appendix II of ref. 10) 
a 


(Aq) 2 = 


Z[q 2 „ - 3(Aq) 2 ] 


( 7 ) 


where a is a nondimensional excitation parameter given by 


a — 


Spp(f) 

S (f ) * - 
PP' 0 L 


: U- -iZ = E^rz=i;3r:=.ii t.:z 


Spp(f) is the spectral density of p at frequency f and S pp (f)i_ is the spectral density of p 
when the r.m.s. value of q equals 1 . Solving equation 7 


( Aq) 2 = 


-f- 


6a/q‘ 


( 8 ) 


« o o q 

Substituting q^ = q + (Aq) and q = 0 from equation 4 into equation 5, the dynamic 
equation becomes 


Aq + 2Q^Aq + Q 2 [3(Aq) 2 - q^]Aq = £2 2 p 


( 9 ) 


From equation 9, the random response is given by (See Appendix II, of ref. 10) 

a ■ 


(Aq 2 ) = 


3(Aq) 2 - q| 


( 10 ) 


Solving equation 10 
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Aq 2 = 


1 + 


Ju^2a/( qj) 


( 11 ) 


The relationship between dynamic response and excitation level is shown in figure 
2. It can be seen that with increasing excitation the rando m response would follow 

4 2 2 

the first solution until it reaches A, where a =%/6, (Aq) = q 0 /6, where it would 

jump to the second solution due to the onset of snap-through motion. The region 
between B and C on the second solution curve is for intermittent snap-through 
motions. Point C, where a = 2qj and (Aq) 2 =q 2 , corresponds to the starting point 

of continuous snap-through motion as found by previous experimental work. It 
should be noted that prediction for points near A and B are subject to uncertainty 
because of the unsteady nature of the intermittent snap-through motion. Figur e 3 

2 

indicates intermittent and continuous snap-through boundaries and values of (Aq) 

2 

for various values of a and q 0 . 

APPLICATION OF SNAP-THROUGH MODEL 

In order to use equations 8 and 1 1 to predict the dynamic response of plates it is 
first necessary to evaluate a, q 0 for different plate parameters and boundary 
conditions using the procedures described below. 

Evaluation of the Non-Dimensional Displacement Parameter 
The nondimensional displacement q = Q/Qn, where Q is the modal 
displacement and Qn is the modal displacement for which the linear restoring force 
due to bending is equal to the nonlinear restoring force due to in-plane stretching. 
The nonlinear effect is significant when Q > 0.3 Qn or q > 0.3. Values of (Q|\|/t) 2 for 
isotropic plates are shown in table 2 for different transverse and in-plane boundaiy 
conditions, as a function of aspect ratio, s, poisson ratio, v, and in-plane stiffness 
parameter, y. It should be noted that (QN/t) 2 is independent of the dimensions and 
elastic modulus of the plates. The steps of derivation for table 2 is found in reference 
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10 and it agrees with other results (ref. 1 and pg. 229 of ref. 5). The significance of y 
is discussed in Appendix I. Based on v = 0.3, a nomograph for finding (OnA) 2 and 
corresponding edge bending strain en for clamped plates is presented in figure 4. 

The Evaluation of the Excitation Parameter 

The non-dimensional excitation parameter a is given by 

s p P < f > 

a=S pp<'c-> L 

Since the non-dimensional force parameter, p, is proportional to the acoustic 
pressure, h, a is related to the acoustic pressure spectral density by 


In reference 10 it can be shown that (see Appendix II) 

S (f ) =64rc 3 c^ fVt 2 Q 2 ki 
hh v o'l o' n 

where c = 0.82 for clamped boundaries and c = 0.38 for simply-supported 
boundaries. 

From the above equations the nomograph in figure 5 has been developed for finding 
a. The use of this nomongraph requires knowledge of the modal damping 
coefficient, C; plate thickness, t; fundamental resonance frequency, f 0 ; plate material 
density, p; acoustic pressure spectral density, Shh(f); and the ratio (Qfs|/t) 2 . An 
approximate value for f 0 can be found for aluminum and titanium alloy plates from 
the nomograph in figure 6. This nomograph is based upon results from reference 5. 
(Q N /t) 2 can be found from figure 4. 

Steps in Using the Design Charts 
The general use of the design charts is outlined in the following steps. 
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1 . Find (QiVt) 2 from figure 4. Then find the bending strain e N corresponding to 
On- 

2. Find f 0 using the figure 6. 

3. Find a from figure 5 using the sound spectral level in dB/Hz. 

4. Find q^ from the compression parameter, R, or from the temperature (ref. 

8). It can also be found by using q 0 = Qo/Qn or £o/en if the static buckling 
displacement Q 0 or bending strain e 0 are known. 

5. The nonlinear response (Aq) 2 corresponding to q^ and a can be found 

from figure 2 or directly calculated from equation 8 and 1 1 . 

6. The r.m.s. strain response is obtained by e r ms = V(Aq) e n- 

7. The values of a corresponding to snap-through boundaries for given 
values of q^ can be found from figure 3. The value of q^ corresponding to 

2 

maximum response (Aq) can be found if the value of a is given. 

The following examples are used to illustrate the steps in using the design charts. 

Example I 

The specimen is an aluminum alloy plate, 12 in. x 15 in. x 0.064 in. with all sides 
clamped and immovable and was tested in the Thermal Acoustic Fatigue Apparatus 
at the NASA Langley Research Center (ref. 16). The modal damping is assumed to 
be 0.02 

From figure 4, using s = 1 .33, (Qi\i/t) 2 is found to be 1 .38 for an immovable edge 
condition. With t = 0.064 in. and a = 12 in., the bending strain on the shorter edge is 
en = 450 microstrain. 

From figure 6, using a = 1 2 in., s = 1 .33, and t = 0.064 in., f is found to be 1 28 
Hz. 

From figure 5, using f = 128, £ = 0.02, (QN/t) 2 = 1-38, t = 0.064 in., Shh = 121 
dB/Hz., the excitation parameter a is found to be 1 .0. 
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The comparison between predicted and experimental results are summarized 
in table 3. For the flat plate, the present method agrees well with a linear method 
(ref. 17) for low SPL but gives a better estimate of the dynamic strain for Shh of 130 
dB'Hz. 

From table 3, the prediction is not very good for the response without snap- 
through for q^O.36. This error is probably due to the presence of imperfections in 

the plate which affects the buckling displacement. However, the predictions for 
continuous snap-through motion are quite good. The effects of small buckling 

displacement becomes less significant compared with the effects of the dynamic 
displacement. For q^=1 .52, the prediction is not very good for intermittent snap- 

through motion. 

Example II 

For this example, the specimen and conditions are similar to Example I except 
that the thickness is 0.040 in. 

From figure 4, (Q N /t) 2 = 1 38. For t = 0.040 in. and a = 1 2 in. the bending strain 
on the shorter edge is en = 195 microstrain. 

From Figure 6, f = 80 Hz. 

From Figure 5, a = 1 when the spectral level is 107 dB/Hz. 

The comparison between predicted and measured results are shown in table 4. 
Agreement was quite good for conditions where qg=0, 0.2, 1.0, 2.48 except for some 
intermittent snap-through. For qQ=2.98, 5.10, 6.7, the snap-through response starts 

at SPL values lower than predicted. Examination of the response spectrum shows 
significant resonance response of the second mode at a frequency close to the 
fundamental frequency. The differences between the results from the present 
method and the finite element method in reference 2 are less than 30 percent. 
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Example III 

For this example the specimen is an aluminum beam, 1 5 in. x 2 in. x 0.032 in., 

clamped at both ends with base excitation of uniform spectral density from 20 - 
100 Hz. The modal damping, £, was found to be 0.017 (details in Appendix III). 
From Figure 4, using s = oo (actually s=6 is adequate) and Sp/Sg = 0.25 (the clamp is 
movable due to bending of the base), it is found that (Qn^) 2 = 1 -7 and en = 71 .8 
microstrain. 

From Figure 6, using s = oo , f = 30 Hz. 

From Figure 5, a=1 when the spectral level is 91 dB/Hz. 

The comparison in table 5 shows very good agreement for the flat beam. The 
large errors for the response of the buckled beam without snap-through are possibly 
due to imperfections. However, the responses with snap-through motion show better 
agreement. 

DISCUSSION 

It is expected that experimental errors included both errors in damping ratio and 
in transverse boundary conditions. The accuracy of the prediction results depends 
mainly on the evaluation of the non-linear stiffness, which is greatly affected by the 
in-plane boundary conditions. 

Near the buckling point, the displacement is greatly affected by imperfection in 
the plate and the two static positions may not be symmetric about the flat position. 
This would have significant influence on the onset of snap-through motions and on 
the non-linear response without snap-through motions. Consequently, the errors in 
response without snap-through were generally greater than errors in response with 
continuous snap-through. It should be noted that intermittent snap-through motions 
are of small magnitude and excited by low excitation level, thus they are generally 
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not as important for fatigue considerations. The continuous snap-through motions 

are more important and can be predicted reasonably well by equation 11. 

For high buckling displacement (q^ > 2.5) there is a possibility of coupling 

between the first two modes in the dynamic response and, therefore, snap-through 
motion may initiate at an excitation level lower than predicted. The (1,1), (2,1) modes 
are found to be coupled during static snap-through (ref. 13). 

For the three examples discussed, the errors of the predicted results are mostly 
within ±30 percent of the experimental results for q^ ^ 2.5 and 35. These ranges 

of buckling and excitation parameters are adequate for most practical purposes. To 
improve the design charts, equation 8 and 1 1 should be derived with consideration 
of higher modes. 

CONCLUDING REMARKS 

A simple and direct design procedure for estimating the nonlinear response of 
flat and buckled plates has been developed. The procedure gives good results for 
moderately buckled plates under moderate excitation. The major source of error for 
highly buckled plates is attributed to the coupling between the first two modes of 
vibration. Other sources of error include imperfection in the plates, damping 
estimates and boundary conditions. For most practical structures, the range of 
application of the design procedure is adequate for first order estimation of the 
nonlinear response. However, the most important step is to evaluate correctly the in- 
plane boundary condition. 
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APPENDIX I 

The Influence of In-Plane Boundary Condition 

From reference 10, the in-plane stresses o x , o y are given by 
Gx = Oxm + Oxc 

Oy = Gym 4- Gyc (A1) 

where G X m. G ym are the mean of g x , G y respectively over the area of the plate, 
o xc , Oy C are the corresponding varying components. 

The strain energy is given by equation A9 of reference 10: 


V = EtabeQ 4 + (a 2 + o 2 - 2vc a ) 

e E xm ym xm ynv 


(A2) 


where e = a constant related to the total strain. energy due to o xc , G yc . 

(see reference 10 or reference 1) 

It should be noted that e depends only on the bending magnitude of the plate and is 
independent of the in-plane boundary condition. G xm , Gym are dependent on the 
boundary conditions and are given by equation A6 of reference 10. 


( 


1 + 


v 


B J 


°xm = 


(1 - v 2 ) 


(Hx + vHy) Q 


1 + 


V 

S B ) 


a = 
ym 


(1 - V 2 ) 


(Hy + vHx) Q 


(A3) 


where Sp, Sp = in-plane stiffness of the plate and boundary respectively. 

H x= 25 b ^ dXdy 


0 0 


Ox 


H = -r^T- J J | dx dy 
y 2ab J 0 J 0 vOy 1 
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The influence of in-plane boundary condition is not considered in reference 1 and 5. 
Note that nonlinear restoring force = KnQ 3 = 5 V e /5Q. 
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APPENDIX II 

Derivation of Reference Sound Spectral Level Spp(f 0 )i_- 

From reference 10 

5C 4C 
S pp(V[_ ~ IT “ jcf 


(B 1 ) 


KQ 


N 


(B2) 


. b . 3 

P=hl where I =J J <J>(x, y) dx dy 
1 1 o o 


(B3) 


K = M il 2 0 = M(2nf 0 ) 2 


(B4) 


M = modal mass = pt I, 


b , a 


= where I = J J <j> (x.y) dx dy 


o o 


Thus K = 4 tc 2 f^ ptl 2 

Substituting P from equation B3, K from equation B6 into equation B2 
h I 


i 


4jt 2 f 2 pt I Q 
o K 2 N 


or h = 4 re 2 f 2 pt Q 


N 


'V 


Therefore, 


W L = 167 r 4 f oP 2t2Q N 




v \) 


S PP^ f o^L 


(B5) 


= 647T 3 cC f 3 0 p 2 t 2 Q 2 n 


(from B1) 


(B 7) 
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c = 


I'J 


APPENDIX III 


Base Excitation Test 

Base excitation involve the excitation of the supporting frame of the specimen 
by electromagnetic shaker. The specimen will experience an equivalent uniform 
pressure excitation proportional to the acceleration of the base and is given by 
(Derivation in ref. 18). 
h = A pt 

where A is the base acceleration in in/sec 2 , p is density of the beam in Ibm/in 3 , t is 
the thickness of the beam in inches. 

Thus Shh(f) = P 2 t 2 S A A(f) 

The spectrum of base acceleration S AA (f) is found to be quite flat above 20 Hz. 

Tests with sinusoidal excitation are performed initially. It is found that the use of 
single mode can predict the fundamental frequency, displacement response 
(measured by a small accelerometer at the center of the beam), strain response (at 
the edge) to within 10 percent of the measured values. 
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Table 3. R.M.S. Strain Response of an Aluminum Plate (12 in. x 15 in. x 

0.064 in.) to Broadband Excitation 

( ) result from reference 17 using linear theory 


Buckling Condition 

Spectral Level 
Shh (dB/Hz) 

Predicted 

erms(microstrain) 

Measured 

erms(microstrain) 

flat q* = 0.0 

100 

35 (30) 

41 

eisi = 450 

110 

115 (95) 

1 1 2 

S L = 121 dB/Hz 

120 

275 (300) 

262 

130 

545 (950) 

443 

buckled q^ = 0.36 

110 

218 st 

85 

e 0 - 270 

115 

270 

173 st 

Ss = 104.2 dB/Hz 

120 

344 

307 

Sc = 115 dB/Hz 

122 

382 

396 

124 

424 

487 


130 

589 

573 

buckled q^ = 1 .52 

105 

41 

31 

e 0 = 553.4 

110 

75 

81 

Ss = 116.8 dB/Hz 

115 

146 

1 33 

Sc = 127.6 dB/Hz 

125 

493 st 

360 st 

130 

620 

556 


st = start of snap-through 




Table 4. R.M.S. Strain Response of an Aluminum Plate (12 in. x 15 in. x 

0.40 in.) to Broadband Excitation 

( ) result from reference 2 using finite element method 


Buckling Condition 

Spectral Level 

Predicted 

Measured 


Shh (dB/Hz) 

erms(microstrain) 

Erms(microstrain) 

a) Flat plate 

110 

159 (190) 

169 

£N — 1 95 

115 

222 

1 84 

S L = 107 dB/Hz 

120 ! 

303 (405) 

265 


125 

410 

382 

buckled q^ = 0.2 

110 

179 (210) st 

252 st 

e 0 =87.2 

115 

236 

326 

Ss = 85.3 dB/Hz 

120 

314 (410) 

414 

Sc = 96.1 dB/Hz 

125 

417 

500 

buckled q^ =1.0 

110 

195 (210) st 

1 60 

e 0 = 195 

118 

291 

206 

Ss = 99.2 dB/Hz 

123 

381 (440) 

470 

S c = 110 dB/Hz 

128 

503 

564 

buckled q^ = 2.48 

105 

78 

89 

e 0 = 307 

110 

225 st 

165 St 

Ss = 107 dB/Hz 

122 

376 

574 

Sc = 117.8 dB/Hz 

127 

488 

638 

buckled q^ = 2.98 

102 

46 

89 

e 0 = 336.6 

106 

78 

165 St 

Ss= 108.7 dB/Hz 

122 

379 st 

470 

S C = 119.5 dB/Hz 

127 

489 

694 

buckled q^ =5.10 

101 

31 

158 

c 0 = 440 

105 

53 

202 st 

Ss =113.2 dB/Hz 

109 

83 

398 

Sc = 124.0 dB/Hz 

119 

356 st 

676 

buckled q^ = 6.7 

101 

26 

78 

e 0 = 505 

105 

42 

140 

Ss = 11 5.8 dB/Hz 

110 

78 

276 st 

S C = 126.6 dB/Hz 

115 

158 

456 


st = i tart of snap-through 



Table 5. R.M.S. Strain Response of an Aluminum Beam (2 in. x 15 in. x 

0.032 in.) to Broadband Excitation 
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Buckling Condition 

Spectral Level 
Shh (dB/Hz) 

Predicted 

erms( mi C r °strain) 

Measured 

erms(nnicrostrain) 

flat 

65 

23.6 

2.4 

£N = 71.8 

76.5 

13.0 

9.2 

S|_ = 91 dB/Hz 

83.4 

25.5 

20.2 

91.7 

49.8 

41 .2 


95.9 

71.1 

61 .0 


102.3 

100.5 

94.0 


103.9 

111.0 

101 .5 


105.5 

122.3 

1 16.0 


108.1 

143.1 

131 .0 


109.8 

158.4 

143.5 

Buckled 

79 

8.1 

2.2 

Po = 3.64 

88.7 

25.7 

7.4 

e 0 = 161.3 

94.1 

60.8 

18.9 

Ss = 94.2 dB/Hz 

95.3 

1 1 1 .9 St 

33.3 St 

S C = 105 dB/Hz 

96.3 

115.2 

46.8 

103.6 

151.4 

1 14.5 


106.4 

172.2 

148.5 


108.0 

186.0 

162.5 


109.4 

199.6 

173.0 


st = start of snap-through 




Force, p 
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w - Q<P(x,y) 


1 R =-1.0 (R u/u c -1) q-Q/Q D 

— 2 R -0.5 M 

3 R= 0.0 

4 R= 0.5 

5 R= 1.0 


Figure 1. Relationship between nondimensional displacement and force. 





(0 N /t) 2 
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Figure 3. Intermittent and continuous snap-through boundaries and relationship 
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t f 

sequence 


Figure 4. Nomograph for evaluating modal displacement, Qn, and surface strain, 
£f\i> when the linear restoring force equals the nonlinear force. 





Random response 



Figure 6. Nomograph for evaluating fundamental frequency, fo- 


continuous 

snap-through 
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